High precision nano scale temperature sensing using single defects in diamond 
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Measuring local temperature with a spatial resolution on the order of a few nanometers has a 
wide range of applications from semiconductor industry over material to life sciences [T]. When 
combined with precision temperature measurement it promises to give excess to small temperature 
changes caused e.g. by chemical reactions or biochemical processes 2 . However, nanoscale 
temperature measurements and precision have excluded each other so far owing to the physical 
processes used for temperature measurement of limited stability of nanoscale probes 3 . Here 
we experimentally demonstrate a novel nanoscale temperature sensing technique based on single 
atomic defects in diamonds. Sensor sizes range from millimeter down to a few tens of nanometers. 
Utilizing the sensitivity of the optically accessible electron spin level structure to temperature 
changes 0] we achieve a temperature noise floor of SmK/v^Hz for single defects in bulk sensors. 
Using doped nanodiamonds as sensors yields temperature measurement with 130mK/v^Hz noise 
floor and accuracies down to 1 mK at length scales of a few ten nanometers. The high sensitivity to 
temperature changes together with excellent spatial resolution combined with outstanding sensor 
stability allows for nanoscale precision temperature determination enough to measure chemical pro- 
cesses of few or single molecules by their reaction heat even in heterogeneous environments like cells. 



Several kinds of nanoscale temperature sensing tech- 
niques have been developed in the recent past [1]. These 
are scanning thermal microscopes (SThM) [5], dispersed 
or scanned individual nanoprobes [3l [6] , direct methods 
like micro-Raman spectroscopy [7] or near-field optical 
temperature measurements [8]. SThMs have tempera- 
ture sensitive elements at a scanning tip (e.g. thermocou- 
ple), the nanoprobes have temperature dependent prop- 
erties (e.g. fluorescence spectrum) which can be accessed 
without direct contact. 

In this study utilize a single quantum system in a solid 
state matrix as a temperature nanoprobe, namely the 
negatively charged nitrogen-vacancy (NV) center in di- 
amond, which allows probe sizes down to ~ 5nm [9|. 
High fidelity control of its ground state electronic and 
nuclear spins has been demonstrated for various quan- 
tum information test experiments [T0"-15 as well as for 
nanometer scale metrology purposes [16.-.19 e.g. measur- 
ing small magnetic and electric fields. Here we show that 
it also allows tracking temperature with high precision. 
Temperature nanoprobes can be either dispersed in the 
specimen to be investigated or used in scanning probe 
geometry (see fig. flk). 

The NV center is a molecular impurity in diamond 
comprising a substitutional nitrogen impurity and an ad- 
jacent carbon vacancy. Optical excitation in a wave- 
length range from 460 nm to 580 nm yields intense flu- 



orescence emission [20 . Excitation also leads to a high 
degree of ground state electron spin polarization (S* = 1, 
the actual sensor level) into its ms = 0(|0)) sublevel 
[2T] . Furthermore the fluorescence decreases upon spin 
flips into the ms = ±1 (|±)) sublevels allowing single 
spin optically detected magnetic resonance (ODMR) ex- 
periments at room temperature [22]. Spins associated 
with the NV center possess long coherence times (~ 1 ms) 
[121 ES] such that coherent control and metrology can be 
achieved with high fidelity. A detectable temperature de- 
pendence of the spin resonance frequencies is observable 
down to temperatures of around 120 K, and very recently 
has been confirmed up to 700 K [24 . 

Optical excitation with 532 nm laser light and fluores- 
cence collection of single NV centers is realized by a con- 
focal microscope which images a bulk diamond sample 
or nanoscale diamonds. Microwave (mw) radiation for 
spin manipulation are applied by appropriate wires close 
to the NV position. The ODMR spectrum of the NV 
center spin usually comprises of two resonance lines cor- 
responding to spin transitions between levels |0) ^ |— ) 
and |0) ^ 1+) [22 (see supporting information). These 
spin levels depend on external parameters like magnetic 
(B) and electric fields, temperature T and strain (see 



fig. fit)) . By analyzing ttie spin Hamiltonian 

H = D{r)Sl + 7NV (B.5, + B^S, + BySy) 

+ Hhf + Hnuc 



(1) 



one can distinguisii between individual contributions. In 
Fig. [ij) the first two largest terms of eq. ^ are illus- 
trated. These are the crystal field which splits spin state 
|0) from states |±) and the Zeeman term due to axial 
magnetic field Bz which splits states |±). The latter 
terms commute. The crystal field parameter D depends 
on temperature T, axial electric field and strain. Un- 
der ambient conditions the temperature dependence is 
cr = dD/dT = -74.2kHz/K [4 . By utilising a ded- 
icated coherent control technique we reduced contribu- 
tions from electric field and strain which are therefore ne- 
glected in eq. (fTl) . The next smaller terms in the Hamilto- 
nian is the transverse magnetic field {B^^ By) which does 
not commute with the previous ones. Transverse fields 
split levels |±) and additionally shift them with respect 
to |0). Therefore they can be confused with temperature 
changes. Sufficiently strong axial magnetic fields sup- 
press small transverse magnetic fields. Eventually, each 
electron spin level is split by hyperfine interaction 



H. 



hf 



E E s.Ai^'i 



(2) 



j=spins i—x^y^z 



to the nitrogen nuclear spin (^^N^ ^ = 1, ^^N^ 
/ = 1/2) [25| and is inhomogeneously broadened by the 
^^C nuclear spin bath [23l [26]. Other minor broaden- 
ing effects are impurity electron spins or fluctuating elec- 
tric or magnetic fields. In diamonds with natural abun- 
dance of ^^C the inhomogeneous ODMR linewidth is 
I/T2* ~ 1 MHz. This has severely limited temperature 
sensitivity in previous measurements where line shifts 
were used to measure temperature, ji) [24l l27l [28] . 

Following the analysis of the spin Hamiltonian we 
adapt our experiments accordingly. Non-commuting 
terms in eq. ([l]), i.e. mainly transverse magnetic fields, 
are suppressed by a sufficiently large magnetic field 
aligned along the NV axis (i.e. B = B^). To suppress 
commuting parasitic effects, e.g. axial electric fields due 
to lattice defects, we use type Ha synthetic HPHT di- 
amonds with low nitrogen content (< 14ppb, see sup- 
porting information). Additionally, we have developed 
a special decoupling sequence (see fig. ^) which elimi- 
nates resonance line broadening due to magnetic inter- 
action with e.g. spin impurities and at the same time is 
sensitive to changes in D (i.e. temperature). 

As illustrated in fig.[l|3 we first create the superposition 
state |0) + |— ) which acquires a phase ipB — ^d during 
time r/2, (e"^^^ |0) + e"^^^ |-)), is then converted to 
g-i^i? |0) — e"^'^^ 1+) and acquires an additional phase 
-^B - ^D during time r/2, (e-^^^^ |0) - |+)). The to- 
tal phase —2if>D = AD • r is proportional to shifts of D 
only. AD is a change of D with respect to an initial 



resonance condition of spin levels and applied mw fre- 
quencies. All intermediate phases (^b caused by quasi 
static fluctuations of magnetic field are canceled. Hence, 
our decoupling sequence is of echo type for magnetic field 
fluctuations and of Ramsey type for changes in crystal 
field D (i.e. D-Ramsey) with a coherence time Td of up 
to ~ 1 ms corresponding to a homogeneous broadening 
of I/Td ^ IkHz. Therefore we increase the phase ac- 
cumulation time r by orders of magnitude from T2* to 
Td which in turn improves the frequency (temperature) 
uncertainty. 

The nuclear spin bath (see eq. Q) does not only lead 
to quasistatic fluctuating magnetic fields but addition- 
ally leads to entanglement between the bath spins and 
the electron spin [26l [29] . The latter effect also appears 
during echo sequences and effectively leads to faster de- 
coherence. The main requirement for this entanglement 
is that Azx and Azy terms of the hyperfine interaction 
are comparable with the nuclear Zeeman energy (see i^hf 
and i^nuc in eq. Q and the supporting information) . We 
face this obstacle by reducing the ^^C concentration for 
the HPHT bulk diamond samples and by increasing the 
magnetic field which reduces the number of bath spins for 
which the above mentioned requirement is fulfilled. The 
effect of these two countermeasures can be estimated by 
numerically simulating the D-Ramsey decay due to the 
spin bath using the cluster expansion method [26 (see 
supporting information). Figure |2]3 shows the expected 
D-Ramsey decay for a few magnetic field and ^^C con- 
centration settings. 

Our experiments are divided into two parts. In the first 
part we test our novel D-Ramsey sequence for single NV 
centers in bulk diamond, to explore the best temperature 
sensitivity achievable. In the second part we apply NV 
centers in nanodiamonds as nanoscale sensors to measure 
temperature changes as a function of distance to a mi- 
crostructure used as local heat source. 

In the first part (see fig. ^) we consecutively mea- 
sured full D-Ramsey oscillations for one day to track 
temperature and possible other effects on D. Figure [2]: 
shows the crystal field parameter D calculated from the 
D-Ramsey oscillations' frequencies. In addition the cur- 
rent temperature was measured by a thermistor close to 
the diamond (see supporting information) and for a value 
cj- = — 78.6±0.5 kHz/K both results coincide (see fig.[2]3). 
The average decay time of the D-Ramsey oscillations is 
T]j = 829±24/is for measurement intervals of 17 minutes 
per point. This corresponds to variations of the crystal 
field parameter D of ±0.2 kHz during one measurement 
interval. Hence, temperature variations during one inter- 
val are at least within ±2.5 mK. Frequency uncertainties 
of fits to the D-Ramsey oscillations yield uncertainties for 
the average temperature during one interval of 1 mK (see 
error bars). For the differences between NV based and 
external sensor based temperature data {S{AD) and ST 
in Scheme [2)3) we obtain a standard deviation of 1 mK. 



Our sensor therefore exhibits sufficient long term stabil- 
ity. 

Furthermore we have investigated the influence of 
laser and microwave heating during this measurement 
sequence by artificially increasing the heating power of 
laser and microwave. Summarizing the microwave radia- 
tion has a negligible effect of much less than 1 mK while 
the laser heats the sample by less than 3 mK on average. 

Figure ^ shows full i^-Ramsey oscillations. Highest 
temperature sensitivity, however, is only achieved for 
longest phase accumulation times r. Taking into account 
only measurements with r = 800 /is we have achieved a 
temperature noise fioor of nj- ^ 5 mK/V^Hz with a max- 
imum measurement rate r ~ 1 Hz. For a given r the 
uncertainty of each measurement point is aj- = n'-f-\fT. 
Noise fioor and the maximum measurement rate depend 
on r approx. as \/l/V and 1/r respectively. 

In a second set of of our experiments utilizes NV 
centers in nanodiamonds. We prepared a glass cover 
slide with a microstructured loop-gap resonator with a 
resonance frequency of ~ 1.8 GHz and a linewidth of 
^ 200 MHz (see fig. pa) which can be excited from a 
distance via radiofrequency (rf) radiation. Afterwards 
nanodiamonds with a nominal size of 50 nm and with 
naturally occuring NV centers (see supporting informa- 
tion) were put in polyvinylalcohol and then spin-coated 
on top of this slide (see fig. Isk). Single NV centers were 
addressed by a confocal microscope and their tempera- 
ture response was investigated. 

Checks of the T^ time and the I^-Ramsey lifetime of 
the nanodiamonds yielded values in the range of 1 to 5 /is. 
Next we heated the sample uniformly to check for the 
homogeneity of the temperature scaling cq- among the 
nanodiamonds. For several tens of NV nanodiamonds 
we found slight variations corresponding to our measure- 
ment accuracy. Next we applied our NV nanoscale ther- 
mometers to confirm heating by the micro structure. To 
this end we measured temperature changes close to the 
microstucture wire using a single NV center in a nan- 
odiamond (see fig. JSJd). When increasing the distance 
from the nanodiamond to the micro resonator by a few 
micrometers the temperature drops until ambient tem- 
perature is reached (see fig. Isfc). 

The temperature noise floor reached in the nanodia- 
mond experiment is nj- ~ 130 mK/V^Hz, i.e. larger than 
for the bulk sample because of the shorter Td value. This 
however allows for faster maximum measurement rates of 
r^ 100 Hz. 

Utilising our novel decoupling scheme nanoscale tem- 
perature sensing method results in a significant increase 
in temperature sensitivities using NV diamond. We prove 
the long term accuracy of the sensor to within 1 mK. In 
addition we demonstrate the applicability of NV nanodi- 
amonds for temperature sensing with similar sensitivities 
recently reported bulk samples [24 . Such nanodiamonds 
can be used as distributed temperature probes to reveal 



spatially resolved temperature profiles with spatial reso- 
lution only limited by their size. 

Temperature sensitivity can be further improved by 
a higher photon collection efficiency, e.g. by solid im- 
mersion lenses [30 or diamond pillar structures ^Slj. In 
addition To can be improved by analogues of higher or- 
der dynamical decoupling sequences like CPMG or UDD 
^32j allowing for micro Kelvin temperature sensitivity. 

Our results pave the way for implementing NV dia- 
mond temperature sensors for material and life sciences. 
One might imagine the integration of NV centers into dia- 
mond AFM cantilever tips [19] . As distributed probes the 
selective attachment of nanodiamonds to specific parts 
of living cells or material and device parts seems feasible 
[33]. While the temperature dependence of D vanishes 
for temperatures below about 120 K other temperature 
dependent properties of the NV center, like the the zero- 
phonon line position of the fiuorescence light, might be- 
come accessible at lower temperatures. 

For temperature measurements with higher time reso- 
lution pump probe measurements can be performed. To 
this end the timing of heat generation and measurement 
sequence need to be properly synchronized. As diamond 
has one of the highest heat conductivities fast response 
is guaranteed. 
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FIG. 1. NV thermometer measurement scheme, (a) Nanodiamonds containing single NV centers can serve as distributed 
probe temperature sensor, (b) Electron spin energy levels are mainly influenced by axial magnetic field Bz and temperature T. 
Bz splits levels |=b) and temperature shifts level |0) with respect to |=b). During free evolution the spin states acquire phases 
(fD, ^B and —LpB- The initial superposition state |0) + |— ) of our decoupling sequence is depicted, (c) D- Ramsey sequence, 
upper part shows mw pulses for spin control and resulting spin states, lower part shows phase accumulation. 
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FIG. 2. NV bulk diamond temperature sensor, (a) Simulation of D-Ramsey coherence decays for ^^C concentration 
c = 0.01 and B^ = 5mT (blue), for c = 0.01 and Bz = 50 mT (orange) and for c = 10"^ and Bz = 50 mT (green). Only 
decoherence due to nuclear spin bath has been considered, (b) Bottom part shows consecutive D-Ramsey oscillations (vertical 
slices) of a single NV center. Horizontal and vertical axis are free evolution time r total time t respectively. The fluorescence 
response is color coded. Upper part shows the crystal field parameter D — 2870.685 MHz + AD deduced from the D-Ramsey 
oscillations (circles) . In addition the current temperature T of the diamond is plotted (yellow curve) . We obtain a temperature 
dependence cr = — 78.6±0.5kHz/K. Middle part shows difference between NV temperature results and additionally measured 
temperature [5 (AD) , 5T). 
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FIG. 3. Nanodiamonds as temperature probes, (a) Schematics of the experiment. Nanodiamonds scattered across a 
glass cover slide. A microstructured loop-gap resonator can be heated by rf radiation. In the region of the current antinode the 
resonator dissipates energy due to ohmic loss, heating the underlying glass slide. Confocal imaging was performed from below 
the glass slide, (b) Lateral confocal microscopy scan showing individual NV centers and part of the loop-gap resonator. The 
local temperature has been measured for three different distances from the resonator (black circles), (insert) Temperature shift 
of a single nanodiamond close to the wire in dependence of the rf frequency for resonator excitation. The resonance frequency 
of the resonator is revealed. 



